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technique together with the harmonic summation and the sector decomposition. 



13th International Workshop on Advanced Computing and Analysis Techniques in Physics Research, 
ACAT2010 
February 22-27, 2010 
Jaipur, India 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonConimercial-Shiu-eAlike Licence. 



http://pos.sissa.it/ 



IR subtraction schemes: integrating the counterterms at NNLO in QCD 



Paolo Bolzoni 



1. Introduction 



In quantum chromodynamics (QCD) and more generally in any quantum field theory with 
massless particles one has to face among others the problem of infrared (IR) divergences when 
computing higher orders corrections. According to the KLN (Kinoshita-Lee-Nauenberg) theorem 
[|l|, §] these IR singularities cancel out once one put together all the contributions at the same order 
which are degenerate to a fixed final state (e.g. a parton state emitting one coUinear and/or soft 
gluon is degenerate to the state without this extra emission). This means that one has to compute the 
sum of several contributions which are separately IR divergent leading to a final physical and finite 
answer. To handle these IR singularities in a general way is non-trivial already at the next-to-leading 
order (NLO) accuracy, where however several solutions are known [|3[, 0, ||, |^ ^ ||]. In recent years 
a lot of effort has been devoted to the extension to the NNLO accuracy [Q [T^, |n|, 12, |l^]. In 
particular the subtraction scheme for QCD jet cross sections defined in |[T^, [T3|] is the one we are 
interested in here. This scheme, initially defined for processes without colored particles in the 
initial state, has been extended to an NNLO-compatible scheme with hadronic initial states [Q]. 
Here we will consider only the case of processes initiated by colorless particles and we briefly 
summarize it as follows. 

In QCD the perturbative expansion for any production rate at NNLO can formally be written 
as 



m jet process we have that the NNLO correction can be 

(1.2) 



dO„^^2-^m+2 + 



d^m+l'^tn+l + 



, VV J 



Considering for example the e+e 
written as 

^NNLO _ 

Jm+2 Jm+l 

where the terms represent the doubly-real, the real-virtual and the doubly-virtual contribution re- 
spectively. Each of them is IR divergent while their sum remains finite. The restriction of the 
phase space to define the physical quantity is realized by the jet functions J„. The basic steps of 



subtraction consist in regularizing all the integrals in Eq.(|1.2D using dimensional regularization in 
d = 4 — 2e; then in reshuffling the singularities between the three terms by adding and subtracting 
suitable counterterms in such a way that we end up with three contributions without IR singularities 
i.e. finite ind = 4 dimensions. In this way Eq.(1.2) becomes 
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In Eq.(|1.4[) above do '2 ' ^^'^ ^^m+2 ^ regularize the singly- and doubly-unresolved limits of 



da^^2 respectively. The last counterterm in Eq.(|L4|) which is da^_^2 '^ must regularize the sinlgy 



RR,A| 



unresolved limits of do '2 ^ ^'^d the doubly-unresolved limits of do '2 ' - Finally in Eq.(l^) we 



have that the counterterms do^j^2 ' ^^id ( /j (iO',„^2 ' ) regularize the singly-unresolved limits of 
do^^^ and /j do^'2' respectively. 

2. Integrating the subtraction counterterms 

In order to complete the definition of the subtraction scheme, one has to compute the integral 



of the counterterms over the factorized one- and two-body phase spaces given in Eqs.(1.5,1.6). 
These integrated counterterms have to be computed md = A — 2e dimensions and the result should 
be given in the form of a Laurent expansion in e. According to the KLN theorem, the e poles of the 
£ expansions of the integrated counterterms have to cancel those in the one-loop correction do^^ 
and the two-loop correction do^^^ giving a finite contribution for the real-virtual and doubly-real 
cross sections. In this proceedings we will discuss the various techniques used to compute the 
Laurent expansions in e for the integrated counterterms. 

According to the different unresolved limits there are different mappings of the external mo- 
menta all of them preserve momentum conservation. The integrated singly-unresolved counterterm 



.RR,Ai 



/i '^'^m+2 ' ^^ Eq.(L5) ( which is also the only one that is needed for a NLO computation) has 



been computed in [15]. As it is shown in [16, 17] the singly-unresolved integrated counterterm 

A, 



RV,Ai , / r j_RR,Ai 



/j [do^jj_^\ ' + ( /i dOj^^^2 ' ) ] m Eq.(L6) can be reduced to three different types of basic integrals: 



these are the collinear, the soft and the soft-collinear integrals . Then the complete counterterm is 



built with them and their non trivial convolutions (see e.g. Sections 2.1 and 2.2 of Ref.[]17|]). As 
far as the doubly-unresolved counterterm S2\^^m+2 ^ ~^'^m+2 '^] ^^ Eq.(L6) is concerned we have 



that the computation of the second one wiU be published in a forthcoming paper [18]. The compu 



tation of the first term is clearly feasible employing the same methods and is a work in progress. 
For J2^^m+2 '^ i*^ turned out that it can be reduced to six different basic integrals divided in cou- 
ples of collinear, soft and coUinear-soft integrals. Analogously to the singly-unresolved case this 
integrated counterterm it is built up with these basic integrals and their non-trivial convolutions. 

To compute the basic integrals that are involved in the definitions of the various integrated 
counterterms we studied different methods. A first fully numerical method for the computation of 
the basic integrals is based on sector decompisition (see [ p^ ] and references therein). To extract the 
poles in e of these integrals a Mathematica package has been written implementing the sector 
decomposition techniques [|^, ^. This program also produces FORTRAN codes directly used by 
numerical integration programs. Another method that has been used to approach the computation 
of these integrals exploited integration-by-parts (IBP) identities. With this method analytical re- 



sults were obtained for some of the singly-unresolved integrals []16[]. However it turned out that the 
method based on Mellin-Barnes (MB) representations is in many cases more efficient and accurate 
than the numerical evaluation via sector decomposition and made the analytic results much more 



feasible than the techniques based on IBP [17]. Both the sector decomposition and the MB tech- 
niques have been implemented independently to cross check all the results (also for those involved 
in the doubly-unresolved counterterm) obtained for the Laurent expansions. 
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To briefly show how the MB method works we want to describe a representative example 
involved in the computation of the /j da^_^2 '^ counterterm. As anticipated a full treatment of it 
will appear elsewhere [[l^]. We begin by considering the following basic coUinear integral 



Ik(x)=x dadv 7^ — -j-T — — -— ] , (2.1) 

where do > 2 is an arbitrary parameter, k = —1,0, 1,2 and a; is a kinematic variable in [0, 1]. Let 



us now consider the particular case do = 2,k = 1. With this choice the integral in Eq.(2.1) is split 
into two different contributions. Of them only the second one produces a pole in e due to the factor 
a^^^^ and we concentrate on it, which explicitly is given by 

/■I a^^^^d — a)"^ 
E(x)=x^ / da-, -, , ,,^ ,^ ' , -^, (2.2) 

where we have neglected an overall factor of r(2 — £)r(l — £)/r(3 — 2e) coming from the inte- 
gration over V. Now applying the well known basic MB formula 

1 1 /■+'" dz 

J^iTW = W)L ^»-'-='=nv+^)r(-z) (2.3, 



two times in Eq.( |2.2| ) we obtain 



E(x)= /+'°° ^^i^z2 y,^^.^e-„„,, r(-zi)r(-z2)r(3-£-z)r(i + £+zi)r(i+z2)r(-£+z 



(27r/)2 r(3-2e)r(l + e) 

(2.4) 



where z = Zi +Z.2- According to the definition of the complex integration path in Eq.(2.3), one 
has to choose the contours in such a way that the poles with a r(. . . +z) dependence are to the 
left of the contour and the poles with a r(. . . — z) dependence are to the right of it. Clearly if 
one starts with contours that do not satisfy this requirement one has to deform them taking into 
account the residua of the integrand while crossing every pole of it. This procedure is automatized 



in the Mathematica package MB . m []21[]. Once this is done this package enables one to easily 
compute the Laurent expansion in e where the coefficients are given as a list of MB integrals. The 
numerical evaluation of the integrals is also implemented in MB . m by use of a simple command 
line. Moreover according to the Cauchy theorem the integration over the complex contour can be 
converted into sums over the residua inside the paths. If we follow these steps for our example 
given in Eq.(2.4) we get 

«=1 \ / m.n=l 



X 



Si(m + n + 2)-Si(m + n)+log(^|) +0{e), (2.5) 

where Si(n) = ^[Lj i denotes the usual harmonic numbers. To obtain an analytic answer the 



harmonic summation involved in Eq.(2.5) has to be performed. In this example, like in many 



other cases encountered in the complete computation of the integrated counterterms, the harmonic 
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summation is feasible and can be performed using the XSummer package [22, E3] in Form 3 . 



W-%. Running a proper script we obtain for our example in Eq.(|2.5|) the following analytic answer: 



1 , ,^, / 1 \ x2(3x2-15x+14) 

E{x) = —+\og{2)[\ ^ ^ ' 



e ^' '\ (l--^)V 2(l-x)2(2-x)2 
(x^ - 9x5 + 33x4 - 78x-n 108x2 - 72x + 16) 

+ - (l-.)3(2-.)3 + "('>■ 0-« 

Note that the limit in x = 1 (corresponding to the collinear limit) is well defined because we have 
that lim_r_>i £" = — 1/e + 53/6 — 161og 2 + 0{e). This concludes the discussion of our example to 
show the mehtod of MB representations. 

As mentioned above, in the forthcoming paper about the counterterm /j dG^^2 '^ we will also 



involve convolutions of the basic integrals, for example the collinear one in Eq.(|2.lD with itself: 



h*Ii{x,y) =y / dadv — — -^ — -^ '— - — ^ '—- 4 (1-a x). 



For this particular integral following the steps described above we find 



5yt _ 1 5/ _ 1 

Ik*Hx,y) = '^^4' 



4,-i5/.-i, , ., (i-4,-i)5/-i , (i-5/_i)4,-i 



2(1+^(1-4,-1)) 2(l + /(l-5/,_i)) 



(2.7) 

1 

+0{e-^). (2.8) 

Here the 0{e^^) contribution is already cumbersome and will not be given here. 

3. Conclusions 



In this proceedings after shortly reviewing the NNLO subtraction scheme developed in [12, 



13, 14] we have discussed the various methods to compute the integrated counterterms focussing 



on the MB representations method. Showing a simple example we described its usage discussing 
some new results: the computation of the integrated counterterm ^2^^^111+2 '^ i^ finished and the 



results wiU be reported in [|18|]. Our methods are clearly applicable to compute the final remaining 
integrated counterterm J2 d'^m+i ^ which is a work in progress. We found that the method based on 
the MB representations can be used to obtain analytic results and in many cases is more efficient 
even numerically than others like sector decomposition. 
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